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ABSTRACT 

Upcoming radio telescopes will allow to study the radio sky at low frequencies with 
unprecedent sensitivity and resolution. New surveys are expected to discover a large 
number of new radio sources, in particular those with a steep radio spectrum. Here 
we investigate the abundance of radio relics, i.e. steep-spectrum diffuse radio emis- 
sion coming from the periphery of galaxy clusters, which are believed to trace shock 
waves induced by cluster mergers. With the advent of comprehensive relic samples 
a framework is needed to analyze statistically the relic abundance. To this end, we 
introduce the probability to find a relic located in a galaxy cluster with given mass 
and redshift allowing us to relate the halo mass function of the Universe with the radio 
relic number counts. Up to date about 45 relics have been reported in the literature 
and we compile the resulting counts, N(> S1.4). In principle, the parameters of the 
distribution could be determined using a sufficiently large relic sample. However, since 
the number of known relics is still small for that purpose we use the MareNostrum 
Universe simulation to determine the relic radio power scaling with cluster mass 
and redshift. Our model is able to reproduce the recently found tentative evidence for 
an increase in the fraction of clusters hosting relics, both with X-ray luminosity and 
redshift, using an X-ray flux limited cluster sample. Moreover, we find that a consid- 
erable fraction of faint relics (5*1.4 S 10 m Jy) reside in clusters with an X-ray flux 
below <3x 10~ 12 ergs -1 cm -2 . Finally, we estimate the number of radio relics which 
await discovery by future low frequency surveys proposed for the Low Frequency Ar- 
ray (LOFAR) and the Westerbork Synthesis Radio Telescope (WSRT). We estimate 
that the WODAN survey proposed for WSRT may discover 900 relics and that the 
LOFAR-Tier 1-120 MHz survey may discover about 2500 relics. However, the actual 
number of newly discovered relics will crucially depend on the existence of sufficiently 
complete galaxy cluster catalogues. 

Key words: cosmology: large-scale structure of the Universe - cosmology: diffuse 
radiation - galaxies: clusters: general - radiation mechanisms: non-thermal - radio 
continuum: general - shock waves - methods: numerical 



1 INTRODUCTION 

Some galaxy clusters show in their outskirts large-scale dif- 
fuse radio emission, which apparently does not originate 
from any individual galaxy. These objects are called 'radio 
relics'. Spectacular examples h ave been found for instance i n 
A 3667 (jRottgering et al.lll997h, A3376 jBagchi et al.|[200r^ . 
and CIZA2242 (|van Weeren et alfcoioh . Diffuse sources are 
difficult to detect due to the low surface brightness and 
due to the steepness of the spectra. Moreover, they can be 
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only classified as relics if galactic foreground and fossil ra- 
dio galaxy emission can be excluded and the hosting galaxy 
cluster is identified. We give a list of currently known relics 
in Section \4. II 



Radio relics show steep spectral slopes, which suggest 
that the origin of the radiation is synchrotron emission. 
Hence, radio relics indicate both the presence of relativistic 
electrons and magnetic fields. There are several approaches 
to estimate the strength of the magnetic field in the region of 
relics. Equipartition arguments have been applied leading to 
field strengths in the range ~ 0.5 — 2 [iG (|Govoni fc Feretta 
2004). For the northwest relic in A 3667 upper limits of the 
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hard X-ray flux in that region pro vide a lower bound fo r 
the magnetic field, namely 1.6 /xG (|Nakazawa et al.l [2009). 
The Rotation Measure distribution of polarized emission 
from sources in the cluster volume or in the background 
allows to constrain the magnetic fiel d strength and spec- 
trum leading to values of ~ 1-5 uG dVogt fc EnfilirJ I2OO3I : 
iKuchar fc Ensslinll2009l ; iBonafede et alj|20ich ~ 

Magnetic fields in galax y clusters are either primordial 
(G rasso fc Rubinstein! l200ll ) or injected in the protoclus- 
ter r egion by AGN and/or galactic winds (|Voik fc Atovanl 
2000). Whatever the origin of the initial seed, some am- 
plification mechanisms are required to account for their 
strength in clusters. Cosmological magnetohydrodynamics 
(MHD) simulations predict magnetic field stre ngths of the 
order of )j,G spread over the cluster volume (e.g. iDolag et al.l 
2005). These studies indicate that the amplification of the 
magnetic field resulting by pure adiabatic contraction is not 
sufficient to explain the observed magnetic field strengths. 
Merger events and accretion of material onto galaxy clusters 
are supposed to drive significant shear-flows and turbulence 
within the intra-cluster medium (ICM). This can in prin- 
ciple amplify ma gnetic fields up to at least /iG levels (see 
IDolag et alj|2008l . and references therein). 

The morphology and temperature distribution of the X- 
ray emission of the clusters which host radio relics indicate 
that relics only occ ur in systems with a n ongoing or recent 
merger, e.g. A 754 (|Macario et al.ll201ll ). For A 3667 it has 
been shown that the relic is located where the bow shock o f 
the moving sub-clump is expected (|Vikhlinin et al.ll200ll ). 
For some other clusters the density and the temperature 
jump of the shock front at the position of the relic have 
been identified (see e.g. iMarkevitchl I2OI0I and references 
therein). This suggest the following scenario for the origin 
of the large-scale radio relics: cluster mergers lead to the 
formation of shock fronts which are responsible of electron 
acceleration causing the relic radio emission. 

Two main mechanisms for the acceleration of elec- 
trons have been proposed to explain radio relics: (i) 
adiaba tic compression of fossil radio p l asma by the shock 
wave (jEnfllin fc Gopal-Krishnal l200ll : lEnfilin fc Briiggenl 
|2002| ) or (ii) diffusi ve sho c k ac ce leration (PSA) by 
the Fer mi-I process (iDrurvl 1 19831 : iBlandford fc Eichlerl 
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Mal kov fc O'C Drurvl |2001). In 
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scenario, ra- 



dio relics should have toroidal and complex filamentary 
morphologies showing very steep, curved radio spectra due 
to inverse Compton (IC) and synchrotron losses. In the 
DSA scenario the electrons are accelerated by multiple 
crossings of the shock front (in a first order Fermi process) 
tracing shocks in the presence of ubiquitous magnetic fields. 
It is worth noting th at other alter native scenarios have also 
been suggested (e.g. Kcshct 2010). 

The formation of radio relics i n cosmologi ca l sim - 
ulations has been studie d e.g. by iHoeft et all i|2008h . 
iBattaglia et~"aH l|2009l ) and ISkillman et alj (|2010l ). The lat- 
ter studied structure formation shocks present in two cosmo- 
logical boxes with a comoving volume of 64 3 h 



" 3 Mpc 3 



and 



20 3 ft" 3 Mpc 3 , and use the non-thermal DSA radio model 
of IHoeft fc Briiggenl (|2007l ) (hereafter HB07) to give an es- 
timation of the radio relic luminosity functi on (RRLF) for 
z = and 1 at 1.4 GHz. On the other hand. ICassano et al.l 
(2010), using a Monte Carlo approach, studied the ocurrence 



of 'radio haloes' in merging galaxy clusters assuming that 
electrons are re-accelerated through MHD turbulence, pos- 
ing interesting constraints for the upcoming Low Frequency 
Array (LOFAR) at 120 MHz. Here we would like to provide 
the appropiate scalings of the relic radiation as a function of 
cluster mass and redshift for a given frequency, as well as to 
give some plausible predictions for upcoming radio surveys 
within the context of the DSA scenario. 

Currently, there are several new radio telescopes un- 
der construction, in particular, for the very low frequency 
regime. Moreover, several existing telescopes are getting sig- 
nificantly improved receivers or backends. For instance, in 
the Netherlands and neighboring countries LOFAR is al- 
most completed. This instrument will survey the sky in the 
frequency range from 30 to 240 MHz. New receivers and 
electronics of the expanded Very Large Array (eVLA) will 
drastically improve the sensitivity of the Very Large Ar- 
ray (VLA). Furthermore, the Westerbork Synthesis Radio 
Telescope (WSRT) will be equipped with focal-plane ar- 
ray receivers which will be optimized for 1.4 GHz observa- 
tions. This will significantly increase the field-of-view and 
will hence improve the survey speed tremendously rising the 
expected number of relic observations. 

In this work we develop a framework to relate the abun- 
dance of galaxy clusters in the Universe to the radio relic 
number counts. To this end we introduce the 'radio relic 
probability density', i.e. the probability of finding a radio 
relic with a given radio power located in a galaxy cluster of 
given mass and redshift. A large sample of observed relics 
would allow to fully determine the probability density func- 
tion for a given halo mass function in the Universe. However, 
since only about 45 relics are presently known we therefore 
use the MareNostrum Universe cosmological simulation 
as a way of determining how the probability density scales 
with cluster mass and redshift. In order to normalize this 
function we compile a list of currently known radio relics. 

From the resulting probability distribution we are also 
able to reproduce the recently found fractions of clusters 
with relics in th e combined NORAS+REF LEX cluster sam- 
ple presented bv lvan Weeren et al.l l|201ld ). Finally, we draw 
some conclusions on the amount of non-identified relics due 
to the fact that the hosting cluster is still not known. More- 
over, we present some estimates on how many relics may 
be identified by upcoming radio surveys assuming plausible 
survey specifications. 

This paper is organized as follows. In Section[2]we derive 
the formalism to estimate the number of radio relics and in- 
troduce the 'radio relic probability function'. We also extend 
the usually assumed sharp transition for the flux detection 
limit to a 'discovery probability' which is more appropriate 
for relic samples. In Section [3] we present the cosmological 
simulation used in this work and briefly summarize the shock 
detection method and the non-thermal radio emission model 
adopted. In Section [3] we present the most up to date ob- 
served relic sample, discuss the normalization of our model 
counts and its comparison to observations and present our 
predictions for upcoming radio surveys. Finally, in Section [5] 
we close the paper with the summary. 

In what follows we assume a flat ACDM cosmology with 
a matter density parameter Q.m = 0.27, an amplitude of 
mass fluctuations as = 0.8 and a Hubble constant Ho = 
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100 h km s" 

Hon]). 



Mpc 1 , with h = 0.7 (e.g. iKomatsu et al.l 



2 HOW TO ESTIMATE THE NUMBER OF 
OBSERVABLE RELICS? 

In this Section we present the formalism aim at estimating 
relic number counts in a given radio survey as well as the 
radio relic probability density. 



2.1 Cumulative radio relic number counts 

How many relics are seen in the sky above a given radio flux 
SV at the observing frequency fobs? The flux of a source 



with luminosity per unit frequency P„ 
redshift z is given by 

S, (py,z) ~ Pv 



AnS(z) 



577(f) located at 



(1) 



where v is the rest frame frequency and di(z) is the lumi- 
nosity distance keeping in mind that the appropiate redshift 
correction for the frequency between the rest and observer 
frames needs to be considered. 

We introduce the luminosity function of 'radio relic clus- 
ters', i.e. the number of galaxy clusters per unit comoving 
volume and logarithmic relic radio power as a function of 
frequency and redshift 



n P (z) := 



d 2 N 



dV c dlogP, 



(2) 



where N is the number of clusters and V c the comoving 
volume. The RRLF is obtained by the convolution of a halo 
mass function (e.g. iTinker et al.l l2008) 



n M (z) := 



d 2 iV 



dVcdlogM 



(M, z) 



with the 'radio relic probability density' of finding a galaxy 
cluster of mass M, redshift z and relic radio power P v . 
Therefore, the RRLF becomes 



n P (z) 



n M {z)p{P v ,M,z)A\ogM 



(3) 



where the relic probability density p(P v ,M,z) fulfills the 
following condition 



p(P„,M,z)d\ogP„ 



(4) 



Hence, integrating Eq. @ allows us to write the total abun- 
dance of relics per logarithmic flux bin as follows 



dN 
d logS 1 ,, 



(&) = 



np{z)——(z)dz. 
az 



(5) 



Note that we have used d log P„ = d log S v since P u depends 
linearly on S v . 

In observations low luminosity radio relics are hard to 
identify since the surface brightness may be too low to ex- 
hibit typical morphological features or spectral index varia- 
tions. Moreover, a diffuse radio object is only identified as 
a relic when the galaxy cluster can be unambigously de- 
tected. Depending on the mass of the cluster this may also 
be challenging. As a consequence, we introduce a 'discovery 



probability', instead of a sharp flux-limit we use a smooth 
transition, which includes both the sensivity of the survey 
and the uncertanties present in the identification. We write 
this probability as 



4>{S V ) = \[ 1 +orf 



logS, - log Sf 



(6) 



where the effective sensitivity, SI , basically gives the flux- 
limit and w the width of the transition. 

Finally, the cumulative radio relic function can be com- 
puted convolving Eq. ((5]) with the 'discovery probability' 
and multiplying by the sky fraction, / s , covered by the ra- 
dio survey 



7V(> logS„)=/ s 



diV 



io g s„ dlogS, 



■(S v )^>(S v )d\ogS v . (7) 



The radio flux-luminosity relation given by Eq. ([T]) and the 
redshift integration of Eq. ([5]) are fully determined by the 
cosmological parameters. Since recent cosmological observa- 
tions show that the resulting parameters are well constrained 
the procedure described above can be considered as a direct 
relation between the radio relic probability density and the 
observed number counts. 



2.2 Radio relic probability density 

If we were to build a perfect radio telescope that could detect 
even the faintest radio emission, which radio power distri- 
bution function linked to structure formation shocks should 
we expect? Merger shocks can persist in the cluster periph- 
ery basically forever, hence, every cluster should show some 
relic radio emission. On the other hand, very bright and 
very faint relics are most likely rare events. We therefore 
expect that there is a typical radio luminosity for a cluster 
with given mass and redshift although the related flux is ev- 
idently below current detection limits. As a consequence, we 
assume that the probability density to find a relic is given 
by a log-normal distribution: 



p{P v ,M, z) oc exp 



(log P v - log P„) 
2er| 



(8) 



where ap is the standard deviation of the logarithmic ra- 
dio power and P v is the mean radio power that scales with 
hosting cluster mass, observed frequency and redshift re- 
spectively. We parametrize this function as follows 



logP, 



log P + Cm x log 



M 



+ C z x log (1 + z) + CV x log 



10 14 - 5 h- 1 M Q 

fobs 



(9) 



1.4 GHz 



The radio power normalization is given by the 'reference 
radio power' Po while the scaling with hosting cluster mass, 
redshift and observing frequency is governed by Cm , C z and 
C v respectively. Formally, a different functional form could 
have been chosen for the radio relic probability density as 
long as the condition given in Eq. Q is fulfilled. However, we 
will show in Section ^. 3l that a log-normal function describes 
reasonably well the radio power distribution of our simulated 
relic samples. 
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3 SIMULATING RADIO RELICS 

In order to simulate radio relics we need to use a galaxy 
cluster sample extracted from a cosmological simulation and 
apply an emission and magnetic model to the present shock 
waves. In this section we present our cosmological simula- 
tion, the method used to detect shocked gas within the sim- 
ulated volume and our radio power emission model. Finally, 
we estimate the parameters of the relic probability density 
using our cosmological simulation. 



3.1 The simulated galaxy cluster sample 

Our simulated galaxy cluster sample was selected from the 
MareNostrum Universe cosmological simulation which 
is a non-radiative hyd rodynamical run of a rep resentative 
region of the Universe l|Gottlober &: Yepes!l2007l ). The sim- 
ulation was run with the smoothed par ticle hydrodynamics 
(SPH) Gadget-2 code dSpringell [20051 '). The adopted cos- 
mology is in agreement with a flat ACDM scenario having 
a matter density parameter S1m = 0.3, a baryon density pa- 
rameter fib = 0.045, an intial matter power spectrum char- 
acterized by a scalar spectral index n — 1 and normalized 
to as — 0.9, and a dimensionless Hubble parameter h = 0.7. 
The simulation started at z = 40 using a linear density field 
represented by 2 x 1024 3 gas and dark matter particles in a 
comoving box of 500 h~ x Mpc on a side. The resulting mass 
resolution for gas and dark matter particles is 8.3 x 10 9 
Mq and 1.5 x 10 9 h" 1 M respectively. 

Identification of bound structures is done using the 
parallel friends-of-f riends (FoF) algorithm described in 
iKlypin et all l|l999l) with a linking length of 0.17 in units 
of the mean interparticle separation. In order to generate 
our galaxy cluster catalogs as a function of cosmic time 
we consider five different redshifts up to z = 1, namely 
2 = 0,0.25,0.5,0.75 and 1, and take the 500 most massive 
galaxy clusters present at each cosmic time. In this way, 
the range of cluster masses we are able to probe goes from 
~ 10 14 h' 1 M Q up to ~ 2.5 x 10 15 h' 1 Mq, meaning that 
the baryonic component of the systems inside the virial ra- 
dius is typically resolved with thousands of gas particles for 
the less massive clusters and with several tens of thousands 
for the most massive ones. 



3.2 Shock finding and radio emission in the 
simulation 

The cosmological SPH code Gadget clearly accounts for 
shock dissipation a s shown by shock tube simulations (e.g. 
ISpringel et aiboOlT ). To this end, artificial viscosity has been 
introduced into SPH, which evaluates the loc al velocity field 
to estimate the dissipation IjMonag han 1992). However, this 
technique is not able to determine the Mach numbers, which 
are needed for combining SPH simulations with paramet- 
ric models for radio emission of relics. Two methods have 
been introduc ed for locating s h ock f ronts and estimating 
their strength: iPfrommer et al.l (I2006T) us es the increase of 
entropy with time while Hoeft et al.l ( 20081 ) evaluates spatial 
entropy gradients in single snapshots of the simulation. Here 
we apply a slightly modified version of the latter approach. 

Briefly, our scheme for locating shock fronts can be sum- 
marized as follows. For a given gas particle we evaluate its 





logPo 


Cm 




C„ 


a P 


This work ('a') 
This work ('b') 


21.35 
21.53 


2.56 
2.22 


3.43 
2.49 


-1.20 
-1.15 


0.85 
0.85 


Skillman et al. 


22.20 


3.65 


3.90 







Table 1. Best fit parameters for the radio relic probability den- 
sity given by Eq. JHJ using our set of MareNostrum clusters for 
magnetic field scaling models 'a' and 'b' (see text). The reference 
radio power, Pq, is obtained in our models using available relic 
observations for normalization (see Section 14.21) . As a compari- 
son shown are the ra dio power scaling parameters obtained by 
ISkillman et al . (2010) in the same redshift range (they assume 
model 'a' and an acceleration efficiency £ c = 0.005). 



pressure gradient and define the shock normal of the parti- 
cle as n = — VP/|VP|. In case that the pressure gradient 
corresponds to a true shock front several conditions must 
be fulfilled. In particular, we demand that (i) the veloc- 
ity field shows a negative divergence, (ii) the density in- 
creases from the upstream to the downstream region and 
(iii) that the latter is also valid for the entropy. Utilizing 
the Rankine-H ugoniot jump conditions for hydrodynamical 
shocks (see e.g. lLandau fc Lifs hitz 1959]) these requirements 
allow to determine the Mach number Mi. For a conserva- 
tive estimate we compute the Mach numbers according to 
all three conditions and then take the minimum. We wish 
to avoid the overestimation of the Mach number since this 
could lead to spurious strong radio emission. We apply this 
shock detection scheme to all gas particles inside a cube of 
size 10 /i -1 Mpc (comoving) centred in the centre of mass 
of the systems available in our FoF catalogs. We consider 
here merger shocks, i.e. shocks introduced by cluster merg- 
ers, which are found to have typical Mach numbers around 
~ 2.5 — 3 (Araya-Melo et al. submitted). We note that also 
fast galaxies in a rather cold ICM may generate shock fronts. 

In order to predict the radio power of the simulated 
shock fronts we need to know the magnetic field strength in 
the downstre am area of the shoc k fronts. Fol lowing our pre- 
vious work in I Hoeft et al.l (|2008l ) and that of Sk illman et al.l 
(2010) we assume that the magnetic field is given by 



B 



B x ( 



10" 



0' 



(10) 



where n B is the local electron density, Bo is a magnetic field 
reference value and r\ is the slope of the density scaling. This 
dependence is motivated by the assumption that in average 
the magnetic field in the ICM is frozen in and that the gas 
motions distribute the magnetic field even to the outskirts 
of the cluster where luminous radio relics are generated. In 
fact, using Faraday ro tation measures in the Coma cluster 
iBonafede etafl (|20ld ) found evidence that magnetic fields 
are spread over the entire ICM. In this work we explore two 
magnetic models. I n the first place w e assume Bo = 0.1 /xG 
and r\ = 2/3 (e.g. iHoeft et al.ll2008f ) which typically leads 
to ~ /J.G values at the ouskirts of galaxy clusters (model 
' a'). W e also adopt the scaling found by Bonafcdc ct al. 
|2010l ) that produces higher magnetic field values (> /iG) at 
these locations (model 'b'). Their best-fit model indicates 
a slightly lower exponent than before but a stronger field 
for an electron density of 10 -4 cm -3 , namely, r\ — 1/2 and 
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log P 14 [Watt Hz" 1 ] log P 14 [Watt Hz" 1 ] log P 14 [Watt Hz" 1 ] 

Figure 1. Radio power probability density for relics inside 3.6 h~ 1 Mpc from the cluster centre at z = and observing frequency 
^obs = 1-4 GHz (magnetic model 'a'; see text). The histograms show results from hosting clusters belonging to the MareNostrum 
simulation for three different mass ranges and the solid lines are log-normal fits according to Eq. JHJ normalized using the condition given 
by Eq. (gjl. 



Bo = 0.8 /iG, respectively. It is worth mentioning that the 
upper limits for the IC emission in the hard X-ray band 
for the northwest relic in A 3667 indicate higher magnetic 
field strength 1.6 /iG; iNakazawa et al.l 120091 ') than ob- 
tained here for a typical electron density of ~ 10~ 4 cm" 3 . 
Hence, this could lead to an overestimation of the hard X- 
ray flux for a similar relic in the simulation. However, it is 
not currently known if the relic in A 3667 hosts an excep- 
tionally strong magnetic field or if relics show in general field 
strengths of a few /xG or more. 

As mentioned in the introduction the emission scenario 
(ii) states that thermal (or mildly relativistic) electrons are 
accelerated at the shock front by DSA. Important evidence 
for this mechan ism comes from the re lic in galaxy clus- 
ter CIZA 2242 (|van Weeren et all [2O10h . In this case the 
observed gradient in the spectral index is consistent with 
electrons accelerated at the shock front, while synchrotron 
and IC losses cause the steeper spectral index in the down- 
stream region. The relic in CIZA 2242 also shows that radio 
emission originates from a rather small volume in the down- 
stream region of the shock with an extent less than 50 kpc. 
In HB07 we have worked out the relation between the radio 
emission and the properties of the shock front and down- 
stream plasma within the context of the DSA model. As- 
suming that the relativistic electron population is advected 
with the downstream plasma and cooled down due to syn- 
chroton looses and IC scattering with CMB photons we are 
able to estimate the total radio emission. In particular, the 
radio power per unit frequency contributed by a SPH gas 
particle i can be written as follows 



P v i = 6.4 x 10 34 ergs" 1 Hz" 1 



Ai 



Mpc 2 10" 4 cm- 3 

s - 3 



0.05 V 1.4 GHz/ V 7keV 

QBd,i//xG) 1+ ^ 
' (BcMB/AtG) 2 + (Bcm/mG) 2 



(11) 



In this formula, Ai represents the surface area given by the 
SPH particle, n e ,i is the electron density, £ e is the electron 
acceleration efficiency, Si is the shock compression factor, 
Td ; i is the post-shock temperature, Ba,i is the post-shock 
magnetic field, Bcmb is the magnetic measure of the CMB 
energy density and ty(A4i) is a function that depends on the 
shock strength. 

We would like to note that the efficiency, £ c , for electron 
acceleration, denotes the fraction of the energy dissipated at 
the shock front that is transferred to supra-thermal parti- 
cles. The lower energy threshold for supra-thermal particles 
is computed from the condition that the power-law distribu- 
tion of supra-thermal electrons must meet the thermal elec- 
tron distribution at the lower energy threshold, see HB07 for 
more details. As a result of this approach the radio emission 
decreases drastically for Mach numbers lower than 3. For 
the computations that follow we simply adopt £ c = 0.005. 
We encourage the reader to see HB07 for details. 



3.3 Estimation of the expected radio power 
scalings using the MareNostrum Universe 

Our aim is to estimate the probability of finding relic ra- 
dio emission coming from an arbitrary galaxy cluster hav- 
ing mass M, located at redshift z and observed at frequency 
fobs- To this aim, we analyse the radio emission produced in 
the different clusters of our synthetic samples. As mentioned 
in Section T3. II we take the 500 most massive clusters at each 
considered redshift (i.e., z — 0,0.25,0.5,0.75,1) and iden- 
tify the shock fronts in each one of them. At all redshifts, 
we evaluate the radio power emitted from cluster relics as 
a function of hosting galaxy cluster mass and observing fre- 
quency. Note that when computing the radio power we con- 
sider all the luminosity caused by structure formation shocks 
within a distance of 3.6 ft." 1 Mpc (comoving) from the cen- 
tre of mass without distinguishing between different relics. 
However, in each cluster, there are typically only one or two 
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prominent relics which contribute most to the final radio 
emission. 

We considered five different frequencies in our analysis, 
namely u obs = 0.12,0.15,0.21,0.325 and 1.4 GHz. Fig. CD 
shows the radio power distribution of relics at z — and 
^obs = 1.4 GHz in the case of magnetic model 'a', where the 
three panels show results for different hosting cluster mass. 
Best fit log-normal functions are also shown. We explore the 
parameter space of our relic cluster samples, given by the 
cluster mass, redshift and observing frequency, and repeat 
the fitting procedure of Fig. [T] (for simplicity we assume 
a constant value for up). In this way, we are able to find 
a set of parameters for the radio relic probability density 
capable to reproduce the mean radio power scalings of our 
synthetic radio relics (see Eq. (JHI). The best-fitting scaling 
parameters for the two magnetic field models adopted are 
shown in Table [3~2l In general, the derived scalings show a 
good agreement. However, magnetic model 'b' displays lower 
radio power scalings with mass and redshift in comparison 
with model 'a', which is most noticeable for the redshift 
evolution. The reason for this can be understood in terms 
of the stronger magnetic field values achieved within the 
context of magnetic model 'b'. Since according to Eq. (|11|) 
the radio emission saturates for large magnetic field values 
the resulting radio power scaling is not so pronounced in 
this case. 

As can be seen from Fig. [T] log- normal functions repro- 
duce reasonably well the radio power distribution of the syn- 
thetic relics. However, since there are possibly more small 
shocks, a better resolution in the simulation may serve to 
alleviate the observed skewness in the radio power distri- 
bution at low cluster masses. Additionally, this could let us 
extend the cluster mass range studied to estimate the mean 
radio power scalings. In the following section we present the 
currently known observed relic sample to further normalize 
our theoretical expectations with observations. 



4 HOW MANY RELICS DO WE EXPECT? 

4.1 Compilation of currently known relics 

By definition, relics are diffuse radio emission in the pe- 
riphery of galaxy clusters without any optical counterpart. 
Hence, relics are commonly searched for by correlating ra- 
dio surveys with large catalogues of galaxy clusters. The 
Westerbork Northern Sky Survey (WENSS) has been car- 
ried out at 325 MHz covering the north sky for declina- 
tions hi gher than 28.5°. The noise level of this survey is 
3.6 mJv (|Rengelink et allll997l ). The NRAO Very Large Ar- 
ray Sky Survey at 1.4 GHz (NVSS) covers the sky north 
of— 40° and has a noise level of 0.45 mjy l|Condon et al.l 
1998). Systematic searches for diffuse radio emission in 
galaxy clusters have been undertaken, for instance, by in- 
specting a sample of 2 05 X-ray bright Abell-t ype clusters 
in the NVSS catalogue (|Giovannini et al1ll999l ). by analyz- 
ing the WENSS data at the position of all Abell clusters 
|Kempner fc Sarazi nl l200ll ). and by searching for steep spec- 
trum sources in th e VLA Low -frequency Sky Survey (V LSS; 
ICohen et aTll2007l ) catalogue jvan Weeren et all [20093 ). 

As described above in more detail, current models for 
the formation of relics are not able to predict the actual 
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Figure 2. Cumulative number of NVSS radio relics. The his- 
togram shows the observed relic sample while the solid and dashed 
lines shows our magnetic models 'a' and 'b' normalized to num- 
ber counts at Sf s 4 = 100 mjy (vertical dashed line). The nor- 
malization leads to logPo = 21.35 and 21.53 respectively (see 
Section l4~2l. 



number of observable relics by themselves, because both, 
the number density of relativistic electrons and the strength 
of magnetic fields, are in general poorly constrained quanti- 
ties. Therefore, we wish to normalize the radio relic number 
counts, N(> S v ), using the number of known radio relics. To 
this end, we have compiled a list of all radio relics reported in 
the literature, as far as we are aware of, which can be seen 
in Table We have included all types of radio relics, i.e. 
Mpc-scale single and double relics in the periphery of clus- 
ters as well as smaller relics inside the cluster volume. A few 
of the small relics might be attributed to the compression 
of fossil radio p lasma (known as 'rad io phoenix' class in the 
terminology of iKempner et alj|2004h . However, we do not 
include phoenixes when computing the relic number counts. 

For each cluster in Table [5] we give the flux of the dif- 
fuse emission which has been classified as 'relic' while the 
contribution of radio haloes has been excluded. In the cases 
where halo and relic radio emission are on top of each other 
due to projection effects we only estimate the flux density of 
the relic emission. In many clusters the diffuse r elic emission 
is fra gmented into multiple pieces, e.g. A 2255 (|Pizzo et al.l 
2008), or shows some prominent p atches and very extende d 
emission as well (e.g. CIZA 2242; Ivan Weeren et aT]|2010h . 
Instead of separating individual relics in a single galaxy clus- 
ter, we combine the flux, S v , of all relics in the cluster which 
is consistent with defining the radio luminosity probability 
for diffuse radio emission in clusters instead of that for relics 
(in the same way as done in Section I3.2|l . For our analysis 
it is not useful to introduce relics as self-contained objects, 
since their identification depends inevitably on observational 
parameters such as sensitivity and resolution. Hence, we give 
in column (3) of Table [2]the entire radio relic flux present in 



© 0000 RAS, MNRAS 000, 000-000 



Radio relic number counts and scalings 7 




Figure 3. Cumulative number of NVSS radio relics for different reshift bins. The histograms show results for the NVSS radio relic 
sample presented in this work, while solid (dashed) lines show the outcome of our magnetic model 'a' ('b'). 



each cluster. To normalize the relic number counts, N(S V ), 
we use the cluster flux at 1.4 GHz because most of the mea- 
surements available are done in the 21 cm band. We also 
estimate the radio power of the relics at 1.4 GHz (see col- 
umn (6) of Table [5} assuming a spectral slope of -1.2 which 
is consistent with the parameter C v given in Table [pT2l Inter- 
estingly, A 3667 displays an outstanding high flux. However, 
this object is not the most radio luminous relic as can be 
seen in Table [5] The ten most luminous relics have fluxes 
5*1.4 > 100 mjy. It is worth noting that several of these lu- 
minous relics have been detected in recent years, namely 
1RXS 06, CIZA 2242, and MACS 0717. On the other hand, 
the faintest relics known to date have a flux of ~ 6 mjy. 

4.2 Normalizing the radio relic number counts 

Basically, we would like to normalize the predicted radio 
relic number counts by using the bright-end of the observed 
number count distribution. As noted above amongst the ten 
brightest relics there are however three relics which have 
been identified only recently. This could indicate that even 
the bright-end of the relics list does not contain all bright 
relic sources on the sky which may introduce an offset in the 
global counts. Hence, in order to estimate number counts 
we assume a fiducial flux of 100 mjy to centre the discovery 
probability. This means that at Sf 4 = 100 mjy half of the 
radio relic emission has been detected. Although this choice 
is arbitrary we take this value as a compromise between 
the lowest and brightest relics in the observed distribution. 
Furthermore, since the lowest flux of known radio relics is 
of a few mjy, we can set w = 0.8, which ensures that the 
discovery probability virtually vanishes below these values. 

The non-detection of a relic could be due to several rea- 
sons. For instance, part of the sky might not be covered by 
deep radio surveys, galactic foreground radiation or bright 
sources in the cluster may obscure the diffuse emission, the 
surface brightness of the diffuse emission may be too low, 
or the related cluster could have not been identified yet. 
All these possibilities for the non-detection of existing dif- 
fuse radio patterns in a galaxy cluster are comprised in the 
complementary discovery probability 1 — <f>. 

Having already a model for the relic discovery prob- 



ability we are now able to normalize the number counts 
to the present observed sample (which we dub as 'NVSS' 
since many candidates have been found by means of that 
survey). In order to do so we take from Table [2] all con- 
firmed relics above a declination of —40° without including 
phoenixes. As can be seen in Fig. the observed number of 
radio relics is well reproduced when using a normalization 
given by logPo = 21.35 and 21.53 for magnetic models 'a' 
and 'b' respectively. However, a degeneration between the 
normalization parameter and the detection threshold exists: 
a higher value for the normalization would imply a threshold 
higher than 100 mjy if one is willing to reproduce observa- 
tions. This would mean that the majority of relics with this 
higher flux has not been detected yet. We consider this pos- 
sibility unlikely. On the other hand, we will show below that 
within the context of the magnetic models considered here 
a lower normalization can be ruled out as a result of the 
analysis of an X-ray flux limited cluster sample and their 
associated relics. The obtained low reference radio power 
(logPo ~ 21.4) is enough to reasonably describe the ob- 
served distributions (see Section T4.3p . Therefore, for the set 
of parameters derived above for the relic radio power proba- 
bility (Cm, C z ,C v ,(jp) we are able to constrain the normal- 
ization very well. 

It is worth noting that to determine radio powers for 
the simulated clusters we had to assume an acceleration ef- 
ficiency, £ c = 0.005 (with £ as defined in HB07 model), and 
we had to assume average scalings for the magnetic field, 
see Eq. ([TU}. For our MareNostrum clusters these assump- 
tions lead to a reference radio power of logPo = 22.23 and 
24.13 in the case of magnetic models 'a' and 'b' respectively. 
Since, as mentiond above, to reproduce observations we re- 
quire lower values for the normalization this implies that the 
acceleration efficiency must be £ c < 0.001. In particular, for 
model 'b' the observed acceleration efficiency could be about 
£ ~ 10" 5 which is more in line with the oretical expectation s 
of DSA in Type la supernova remnants (jEdmon et al ] |2011f l. 
Further increase of the magnetic field values, as suggested by 
iNakazawa et alj (|2009t ). would reduce the required efficiency 
even more. 

Most of the observed radio relics have a redshift lower 
than 0.3. In the sample there are only five relics with higher 
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Figure 4. Estimated relic radio power at 1.4 GHz as a function of X-ray cluster luminosity (left panel) and temperature (right panel) 
for the NVSS relic sample. All X-ray quantities are derived in the ROSAT 0.1 — 2.4 keV band (see Tabled. 



redshift, and only one of them is located at z > 0.5. We 
wish to compare these numbers to the predictions according 
to the radio relic probability distribution. We simply split 
the result into the redshift intervals given above, as can be 
seen in Fig. [3] Apparently, our models predict more relics 
than observed for z > 0.3. This might indicate that relics 
in distant clusters are more difficult to detect (e.g. due to 
resolution effects), that clusters would need to be more X- 
ray luminous to be found, or that our scaling parameter C z 
derived from the simulations does not agree with the actual 
redshift evolution. In particular, model 'b' seems to better 
reproduce observations at all redshifts. However, it is im- 
portant to realize that the number of both predicted and 
observed relics with z > 0.3 is very small, so we should be 
cautious with any interpretation. Much more extensive cat- 
alogues of relics are needed to draw a significant conclusion 
about the redshift evolution. 

The most distant cluster which hosts a relic is 
MACS 0717. Fig. [3] indicates that in our model the high- 
est redshift relics should have fluxes Si. 4 within the range 
10 - 50 mjy. Instead, the relic in MACS 0717 has a flux 
of about 140 mjy, indicating that this system is an out- 
standing radio relic. In fact, it is the most luminous relic 
known to date with P1.4 ~ 2 x 10 26 W Hz -1 . For instance, 
using scalings resulting from magnetic model 'a', the mean 
relic radio power of clusters having the mass and redshift 
of MACS 0717 is P1.4 S 4.1 x 10 24 W Hz" 1 (see Eq. ©). 
Hence, the luminosity of the relic is about 2a p higher than 
the mean relic luminosity, so it is a rare event but still rea- 
sonably likely considering all clusters in the Universe. 

In general, we can study the deviation between the 
estimated radio power from observations and the ex- 
pected mean radio power at a given redshift and cluster 
mass. We can quantify this deviation in terms of Ap := 
log (P1.4/P1.4) /up, which measures the difference between 
the logarithmic radio power estimated from observations 
and the peak of the radio relic probability distribution in 
units of the parameter a p. In order to make a simple esti- 



mate we adopt the cluster X-ray luminosities given in col- 
umn (8) of Table [2] to compute the cluster masses. In what 
follows we use a L x — M relation similar to that given by 
iPratt et al] (|2009l ) that will be presented in Section 031 In 
column (7) of Table [2] we give Ap for all relics in our ob- 
served sample adopting the scalings derived with the mag- 
netic model 'a' to estimate the mean radio power. As ex- 
pected most of the relics display a mean deviation of ~ 2<rp 
since we are observing the brightest (or close by) relics in 
the sky. Some of the relics show an unexpected high devia- 
tion from the mean radio power, Pi. 4, given by the hosting 
cluster mass and redshift. This may serve as an indication of 
the need to further investigate these systems in more detail 
to confirm their relic nature. However, the large deviations 
may also come from uncertanties in the cluster mass esti- 
mate. We do not pretend here to give a rigorous derivation 
of the hosting cluster masses but only assess the global ten- 
dency of the sample. 



4.3 The X-ray — radio power relation 

For giant radio haloes in massive galaxy clusters a close 
correlation between r adio power and X -ray luminosity 



has been found (e.g . Venturi et al.l l2007t ). In particular 



lEnfilin fc Rottgerind ||2002T ) suggested that the radio power 
of observed haloes scales with cluster X-ray luminosity ac- 
cording to 



10 24 W Hz~ 



10 45 erg s" 1 



(12) 



with parameters a v — 5.36 and b v — 1.69. In contrast, the 
radio power of relics show a large scatter for a given X-ray 
luminosity or temperature as can be seen in Fig. [4] This fact 
precisely reflects our starting point, namely, the recognition 
that the radio power of relics varies strongly for a given 
galaxy cluster mass. Therefore, this motivated us to intro- 
duce the radio power probability distribution. Formally, we 
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Figure 5. Luminosity function of radio relics for models 'a' 
and 'b' (solid and dashed lines respectively) and radio haloes 
(dot-dashed line) at 1.4 GHz and 2 = 0. The radio halo 
luminosity function i s an analytic approximation taken from 
lEnfilin fc Rottgerinel fe002l 1 under the assumption that a con- 
stant fraction f r ^ = 1/3 of the clusters contain radio haloes. 




9 



S x [erg s 1 cm" 2 ] 

Figure 6. Cumulative fraction of clusters which host diffuse relic 
emission. The cumulative number is depicted as a function of X- 
ray flux, measured in the ROSAT 0.1 — 2.4 keV band. The curves 
give the cumulative number for a relic flux of 1, 10, and 100 mjy. 
In addition the completeness limit of the REFLEX cluster sample 
is given. 



can relate the mean radio power, Pi. 4, to the X-ray luminos- 
ity of clusters at z — using the Lx — M relation in Eq. ([9]) 
for the parameters given in Table I3T21 Comparing this result 
with Eq. (|12[) we get for a v and b v the values 1.3 x 10~ 3 
and 1.5, respectively. Interestingly, we find a similar expo- 
nent but a much lower proportionality constant. This seems 
to indicate that there are much less bright radio relics than 
haloes. However, we have to keep in mind that when com- 
puting the RRLF not only the mean radio power but the 
radio power distribution function needs to be taken into ac- 
count, see Eq. (|3]). As a consequence we expect more radio 
relics than haloes as can be seen in Fig. [S] 



4.4 Radio relics and an X-ray selected cluster 
sample 

The radio power probability density introduced above allows 
us to predict the fraction of clusters in an X-ray selected 
sample which host a radio relic. As a first step, we introduce 
the differential distribution with respect to cluster X-ray flux 



d 2 N 



d log Sx d log S v 
In a similar way as done in Eq. ([5]) we can write 



(13) 



/is Y s = / n M (z) p(P„, M, z) dlo * Sx dz. (14) 

Note that radio power is a function of radio flux and redshift, 
similarly, the X-ray luminosity is a function of X-ray flux and 
redshift. Using the Lx ~ M relation (see below) to estimate 
the cluster mass, based on its X-ray luminosity, we can write 
p{P u ,M,z) = p(P v (S v ,z),M(Sx,z),z). We now introduce 
an X-ray flux threshold, Sx an d assume that only clusters 
with a flux above the threshold are detected. This allows us 



to determine the cumulative fraction of clusters with Sx > 
Sx that host diffuse relic emission with a given flux S v . 
Integrating the previous equation leads to 



N X {S V ) 



nsxSv dlogSx, 



(15) 



where the normalization factor, Nx{S„), is determined by 
the condition Fx(> 0,S U ) = 1. Fig. HJ] shows the cumu- 
lative fraction of clusters for three different radio fluxes 
(for the sake of simplicity we assume only the radio power 
scaling parameters that result from model 'a' throughout 
this section; adopting those of model 'b' do not modify our 
main conclusions). About 80% of the clusters which host 
diffuse relics with 100 mjy have an X-ray flux larger than 
3 x 10~ 12 ergs -1 Hz -1 , which corresponds to the complete- 
ness limit of the REFLEX cluster sample (|Bohringer et al.l 
2004). Since candidate radio relics are commonly identified 
by cross-correlating radio and X-ray catalogues, the NORAS 
and REFLEX cluster catalogues are well suited to identify 
luminous relics. In contrast, for faint relics of about 10 mjy 
only ~ 40% of the hosting clusters are expected to have 
fluxes above the REFLEX X-ray flux limit. This means that, 
if upcoming surveys will allow the detection of diffuse radio 
structures with fluxes of about 1 mjy, significantly deeper 
X-ray cluster catalogues will be needed to identify the ma- 
jority of radio relics. 

In a similar way we can determine the cumulative frac- 
tion of clusters which host diffuse relic emission as a function 
of radio flux. To this end we introduce the fraction of galaxy 
clusters with detectable relics at a given mass and redshift 



U(M,z) 



P {P„,M,z)<j ) {S v )d\ogP l ,, (16) 



where <f> is the discovery probability introduced earlier. Note 
that for w — v we can mimic a Heaviside-function, i.e. only 
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relics above SI are detected. We can now determine the 
cumulative fraction of clusters hosting radio relics per X- 
ray luminosity bin, A log Lx , as follows 

1 dlogM dV c 



F v {> ST,L X ) 



N V {L X ) J nMJ 'dlogL x dz 

where nuj '■= riM f<t>{M,z) and the normalization fac- 
tor, N v (Lx), is given by the condition F v {> 0, Lx) = 1. 
Fig. [7] shows the cumulative fraction for different cluster X- 
ray luminosities. As expected, only a small fraction of clus- 
ters show radio relics for current detection limits of about 
10 mjy. The fraction decreases strongly with the cluster X- 
ray luminosity. For instance, ~ 20% of clusters with an X-ray 
luminosity of about 3 x 10 45 ergs -1 are expected to host a 
relic with a flux of 10 mjy or brighter, while only ~ 0.3% 
of clusters with 3 x 10 44 ergs -1 are expected to do so. Note 
that in Eq. (|17l) we have assumed that all clusters with given 
X-ray luminosity can be detected. To compare the result to 
X-ray selected cluster samples we need to introduce an X-ray 
flux limit. As a result, a large number of faint relics residing 
in distant clusters fa lls below the flux l i mit. 

In a recent work Ivan Weeren et al. I J2011ch selected 544 
clusters fr om the NORAS ( Bohringer et all 1200(1 ) and the 



dz (17) 



REFLEX (Bo hringer et al.l 20041 ) cluster samples with an 



X-ray flux above 3 x 10 - ergs - cm - and located outside 
the galactic plane. Up to this flux the REFLEX sample is 
virtually complete. On the other hand the NORAS sample 
is almost 50% complete. Interestingly, these authors show 
that 16 out of the 544 clusters of the combined list contain 
at least one radio relic and found evidence for an increase of 
the fraction of clusters which host relics with cluster X-ray 
luminosity and redshift. 

Eq. (JTTJ) allows us to determine the fraction of clusters 
with relics. The X-ray flux limit imposes an upper limit in 
the redshift integral, i.e. allowed redshifts must fulfill the 
condition z < z(S t x , L x ). First we wish to reproduce the 
dN c \/dL x and dN c \/dz dist r ibution s of the cluster sample 
selected in lvan Weeren et all (|201ld l. To this end we rewrite 
Eq. (JTTJ) as follows 

dlogM dVc 

dz 



dJVci 
dL x 



riM 



dL> 



dz 



(18) 



and 



diVci 
dz 



n M 



dlogM dV c 



dLj 



dz 



dL x 



(19) 



where N c \ is the number of clusters and / s indicates the sky 
fraction covered by the selected cluster sample, we estimate 
it to be 359^0]. The integration boundaries, zth(<Sjc > Lx) and 
LxiSxjZ) axe obtained from the flux limit in the survey, 
see Eq. |T]). To perform the integration we need to relate 
cluster mass and X-ray luminosity. We assess the cluster 
luminosity and redshift distributi ons by choos i ng an appro- 
priate Lx (M) relation. Recently, IPratt et al.l (|2009T l inves- 
tigated cluster scaling relations in detail. They provide best 



1 The sky fraction is estimated by f^y X /; nt X (1 — / ga i), where 
/sky is the sky fraction covered by the survey which overlaps 
NVSS, /i nt is the completeness of the survey, and / ga j is the 
fraction of clusters located at a galactic latitude lower than 20° . 
For NORAS and REFLEX we assume 50%, 50% and 33%, and 
34%, 90%, and 33%, for these quantities respectively. The sum of 
the two contributions leads to the estimated value. 
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Figure 7. Cumulative fraction of clusters with relics as a func- 
tion of radio flux. The cumulative number is depicted for four 
different X-ray luminosities. 



fit parameters for Lx.soo — M500, where both quantities are 
measured within J?5oo- However, our sample differs in sev- 
eral respects to theirs: (i) the halo mass function, um, is 
given for M200, (ii) we do not extrapolate the luminosi- 
ties to -R500 , (iii) our sample extends up to z ~ 0.5 while 
IPratt et al.l (|2009l ) use only clusters with z < 0.2, and (iv) 
we use a slightly different value for Qm- Adopting the scaling 
relation 



5 x lO^ergs- 1 



= E(z) 



7/3 



M 2 , 



3.4 x 10 14 h- 1 M© 



(20) 



with E{z) = ^/Qm(1 + z) 3 + Q A we are able to reproduce 
the observed distributions reasonably well (see Fig. [8]). 

To compute the abundance of clusters which actually 
hosts a detectable relic we have to use umj = «m/(M, z), 
instead of um in Eqs. (fl8|) and (fi"9)l . Fig. [8] (red solid lines) 
shows the resulting X-ray luminosity and redshift distri- 
butions. Note we use here a lower effective sensitivity and 
smaller width in <f>(S v ) because the clusters are known. As 
discussed above we may miss the discovery of a bright relic, 
since the cluster has not yet been identified. For instance, 
the double relic in PLCK G287.0 has been discovered only 
after the detection of the hosting cluster with Planck satel- 
lite even if the diffuse emission is clearly visible in NVSS. 
Since we consider here the relics in the NORAS+REFLEX 
sample all clusters are known by construction. We model 
the discovery probability in the following way: (i) we argue 
that the width is smaller than for the overall sample, we 
take w = 0.2, (ii) we adjust SfH to reproduce the fraction of 
clusters with relics found in the sample, using Sf 4 = 27mJy 
we obtain a fraction of 3%. Hence, the effective sensitivity 
adopted corresponds to 60 times the r.m.s. noise level in the 
NVSS survey. In this way we find that the fraction of relics 
in an X-ray flux limited cluster sample should indeed in- 
crease with X-ray luminosity and redshift as shown in Fig.E 
The redshift distribution is a crucial test for the mean radio 
power scaling parameter C z - 



© 0000 RAS, MNRAS 000, 000-000 



Radio relic number counts and scalings 11 




Figure 8. Luminosity (left panel) and redshift (right panel) distributions of X-ray clusters in the NORAS+REFLEX sample with and 
without radio relics (red and blue histrograms respectively). Blue and red solid lines show the distributions expected in our model while 
the black solid line shows the fraction of clusters with relics (ratio between the red and blue solid lines). The green dashed lines show 
the expected fraction of clusters with relics for the NORAS sample in the upcoming LOFAR-Tier 1-120 MHz (see Section 14.511 . 



4.5 Predictions for upcoming surveys 

In Table [3] we have summarized specifications for planned 
surveys with LOFAR and the WSRT. We would like to 
give some plausible estimates for the expected number of 
relics to be discovered by upcoming radio surveys. However, 
we have to remember that there may be several uncertain- 
ties difficult to quantify. As we already mentioned the de- 
termination of p(P v ,M,z) is affected by the limitations of 
our simulation and by the adopted physical model used to 
relate the Mach number with the relic radio power. How- 
ever, since we are able to reproduce the trends found for 
the NORAS+REFLEX sample we conclude that our ap- 
proach has resulted in a reasonable set of parameters. In 
the present paper we use the radio relic probability density 
estimated from the MareNostrum simulation and leave for 
future work a more extensive modelling of the radio power 
emission. 

To compute number counts for future surveys we also 
need to assess the discovery probability for each survey. 
Since it is beyond the scope of this work to model this 
in detail we simply adopt a conservative approach: we as- 
sume that the NVSS detection parameters hold similarly for 
the upcoming surveys, i.e. we take w — 0.8 (detection/non- 
detection transition of the instrument) and b := S„ ja v ~ 
200 (ratio of the effective sensitivity of the overall relic sam- 
ple to the survey noise) as our fiducial parameters. The last 
condition means that for half of the clusters which host dif- 
fuse radio emission with a flux above 200 times the noise 
level of the survey a radio relic will be detected. At this 
point, an important remark needs to be done in relation to 
the sensitivity per beam achieved by a given radiotelescope. 
The next generation of radio surveys will presumably in- 
crease their beam resolution at least by a factor of a few. 



In principle, for some of the relics, this could imply the re- 
quirement of higher b values than assumed here which may 
lead to an overestimation of the predicted number counts. 
For instance, in its final configuration the LOFAR telescope 
is expected to reach an angular resolution of ~ 5 arcsec. 
However, since this instrument has many short baselines it is 
always possible to smooth images down to typical NVSS res- 
olution values (i.e. ~ 45 arcsec) without increasing too much 
the resulting r.m.s. sensitivity. Nevertheless, we have to keep 
in mind that the final predicted number counts will depend 
on the adopted radio power scalings and detection parame- 
ters. For instance, in the case of the LOFAR-Tier 1-120 MHz 
survey, if we let the effective detection threshold to vary in 
the range b — 150 — 300, keeping the remaining parameters 
fixed, our predictions will increase (decrease) by a factor of 
two (a half) with respect to the prediction corresponding to 
b ~ 200. Similarly, one could assess the impact of varying 
some of the radio power scaling parameters keeping the rest 
unchanged. In particular, if we let the slope of the scaling 
with cluster mass to vary in the range Cm = 1-5 — 3.5 we 
get brighter (fainter) relics located in clusters with masses 
below (above) 10 14 ' 5 ft -1 Mq. The predictions in this case 
will also increase (decrease) in a similar amount as before. 

For calculation purposes here we use the radio power 
scalings derived from the simulation assuming the fiducial 
detection paramaters presented above. Table [3] gives the to- 
tal number of expected relics up to z = 0.3, as well as relics 
with 0.3 < z < 0.5, 0.5 < z < 1, and z > 1 for mag- 
netic models 'a' and 'b' (in brackets). For z > 1 models 
with higher magnetic fields would generally produce less ra- 
dio relics as a consequence of the resulting redshift evolu- 
tion. Under these assumptions we expect that the LOFAR- 
Tier 1-120 MHz survey and WODAN large sky coverage 
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survey should reveal several thousands of radio relics, due 
to the huge improvement in sensitivity that will presum- 
ably be achieved by these surveys. However, a given survey 
may provide candidates for radio relics, by cross-correlation 
with positions of known galaxy clusters. This means that 
deep follow up observations may be needed to confirm ra- 
dio relics in the clusters. Interestingly, within the context of 
our simple scaling for the magnetic field, we also expect a 
significant number of relics with z > 0.5 and z > 1 to be 
detected. The actual number of relics at high redshift will in 
particular serve to constrain the redshift evolution of mag- 
netic fields in clusters allowing us to further refine our model 
prescriptions. 

As noted above, the number of relics detected in up- 
coming surveys will crucially depend on the cluster database 
used to correlate candidates displaying diffuse radio emission 
with known cluster positions. As a consequence, the final 
number of unambigously identified relics could be below our 
model expectations. To quantify this we estimate how many 
relics LOFAR should find in the Tier 1-120 MHz configura- 
tion for the NORAS cluster sample introduced before. Note 
that the LOFAR survey and the NORAS sample cover both 
the north sky. Based on the effective sensitivity found in the 
previous section for relics in this galaxy cluster sample we 
adopt here a flux threshold 60 times that of the r.m.s. noise 
level of the LOFAR-Tier 1-120 MHz survey {Sf 4 /a T i CI i), 
which results in an effective sensitivity of 6 mjy (see dashed 
lines in Fig. [Sj . We find that LOFAR should discover relics 
in more than 50% of the clusters. In the case of luminous 
clusters the fraction can be as high as 90%. 

Using Eq. (|15|l we can also estimate what is the required 
sensitivity in X-ray surveys to find at least a fraction of the 
relics that can be potentially discovered in a radio survey. 
Fig. [9] indicates that for relics with 20 mjy in the LOFAR- 
Tier 1-120 MHz survey, about 50% of the relics might be 
identified by cross-correlation if the X-ray surveys are com- 
plete up to at least 4 x 10~ 13 ergs~ 1 cm~ 2 (in doing this 
calculation we have used magnetic model 'a' as done in the 
previous section). Hence, an all-sky X-ray catalogue with 
an X-ray flux limit one order of magnitude below that of 
the REFLEX sample is necessary to identify a considerable 
fraction of the relics. 



5 SUMMARY 

Radio relics are believed to trace merger shock fronts in 
galaxy clusters. The radio luminosity of shock fronts de- 
pends strongly on the Mach number of the shock, but also 
on the size of the front and on the magnetic field present 
in the downstream region. Even if in every cluster there are 
shock fronts related with past merger events, the actual ra- 
dio luminosity caused by the shocks may vary strongly, from 
no detectable radio emission at all to the presence of lumi- 
nous radio relics. To describe the large spread of radio lumi- 
nosities more formally, we have introduced the radio power 
probability distribution, p(P v , M, z), aim at assesing the like- 
lihood of relics in galaxy clusters with a given mass and 
redshift in a given frequency. We use the MareNostrum 
Universe simulation to estimate the probability distribu- 
tion. To this aim, we selected the 500 most massive clusters 
at 5 different redshifts up to z = 1 to detect shock fronts 




S x [erg s 1 cm" 2 ] 

Figure 9. Same as Fig.[6]but with specifications for the LOFAR- 
Tier 1-120 MHz survey: v \, s = 120 MHz and an effective sensi- 
tivity S^ff = 20 mjy. The curves give the cumulative number for 
a relic flux of 2, 20 and 200 mjy. As a vertical dashed line is 
also shown the X-ray completeness limit of the REFLEX cluster 
sample, Sr, = 3 X 10~ 12 erg s _1 cm~ 2 . 

in assembling galaxy clus t ers. T hen, we apply the scheme 
developed by iHoeft et alj (|2008T ) for estimating their radio 
relic luminosity. Based on the distribution of radio relic lu- 
minosities of the simulated clusters we conclude that the ra- 
dio power probability is well approximated by a log-normal 
distribution. Moreover, using our galaxy cluster samples, we 
are able to estimate how the radio relic distributions scale 
with cluster mass, redshift, and observing frequency. 

Using the radio power probability distribution we wish 
to determine the relic number counts, N(> Su)- Basically 
this is given by a convolution of the probability distribu- 
tion and the dark matter halo mass function. However, radio 
relics are not straightforwardly identified in radio observa- 
tions, therefore, even luminous relics are possibly present 
in radio catalogues, but not yet identified as relics. For in- 
stance, the relics in 1RXS 06+42 and in CIZA 2242 are 
bright systems present in the WENSS catalogue, but have 
only recently been reported as relics. We therefore intro- 
duce the discovery probability as a function of radio flux, 
4>{S V ). The number counts are obtained by a convolution of 
all three, the halo mass function, the radio power probability 
distribution, and the discovery probability. 

It is important to remark that it is not possible to pre- 
dict radio flux number counts of relics purely from cosmolog- 
ical simulations. In this regard a major source of uncertainty 
is the efficiency of the electron acceleration at the shock 
front. We therefore use the observed relic number counts to 
determine the reference normalization for the radio power 
probability distribution. 

The resulting framework allows us to estimate the num- 
ber of detectable relics in upcoming radio surveys. In the 
following we summarize our main conclusions: 

• To evaluate the MareNostrum Universe simulation 
we assumed an electron acceleration efficiency, £ c = 0.005 
and two different magnetic field scalings with local electron 
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(2008) 



G91: lGiovannini et al] 199j ). G09: IGiovannini et al 
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12009 
l201lh . 



G10: iGiovannini et all ll2010l) , 
K07:lKocevski et alj j2007t). 



Giacintuccietal^ 

GF00: IGiovannini fc Ferettil ll2000l) , GoOl : iGovoni et alj d200ll), G oll: lGovoni et alj 

KS 01: lKempner fc Sarazinl ll200ll ), M9 8: iMarkevitch et alj jl998l) , M03: iMarkevitch et al 1 1|2003|). Mall: iMacario et alj 
201ll).Op: Og r ean e t al. (in p rep.), OlPlOgrean et alj ll201ll) . P0 4:|Popesso et alj ll2004h, P D09: |Pizzo fc de Bruvnl (120091) 



V07: IVenturi et al 



greanet al. (in prep.j, Ull: Ugrcan et al. (2011j, .FU4: r'opcsso et alj I2UU4I). F 1J09: Fizzo & dc ri ruvnl 1121 
.1 1*20071). yW09 a: Ivan Weeren et alj '2009a[), vW09 b: Ivan Weeren et alj '2009b|) vW09 c: Ivan Weeren et al 



'2011ch 



. 3093) ■ 

vWl O: Ivan Weeren et all <201(fl . vWlla: lvan Weeren et alj <201lj). vWllb: lvan Weeren et alj ll2011rj), vW llc: lvan Weeren et alj 
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vWp: van Weeren et a l. (in prep. ), RB02 : lReiprich fc Bohringerl j2002t), R97 : IRottgering et alj U997I). R 99: lReid et alj '199&I) . 
S01: ISlee et"al] '200ll ). Sh0 8: IShen et alj '20081). Su03 : ISubrahmanvan et alj 1I200I ). V08: IVenturi et alj J200S|) . Vll: IVenturi et alj 
j201ll ). ^Giovannini et al] j201ll ) has classified this source as a radio halo. 

Table 2. List of currently known relics extracted from the literature. Columns: (1) Cluster name, (2) Classification (R: single relic; D: 
double relic; P: phoenix; C: diffuse radio emission detected [more observations are needed to confirm its nature]; F: probably misclassified 
as relic), (3) Radio flux (all relics in the cluster are considered), (4) Observed frequency, (5) Redshift, (6) Radio power (computed 
assuming a spectral index of —1.2), (7) Radio power deviation (model 'a'; see text), (8) Cluster X-ray luminosity, (9) Cluster X-ray 
temperature, (10) Checkmark if within NVSS relic sample (8 > -40°), (11) Checkmark if within NORAS+REFLEX cluster sample (in 



brackets if below flux limit Sr = 3 X 10 
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Table 3. Properties of various present and upcoming radio surveys. For LOFAR and WODAN different configurations are also shown. 
The colums correspond to: name of the radio survey, observing frequency, survey noise level, observed sky area (or declination limit), 
corresponding sky fraction and approximate number of expected relics using our set of fiducial parameters as a function of redshift for 
magnetic models 'a' and 'b' (in brackets). It is worth noting that these are only plausible estimates for upcoming radio surveys but are 
not meant to be definitive values (see Section l4.5H . 



density in the simulation. Normalizing the radio power 
probability distribution by the list of known NVSS relics 
resulted in lower values for the reference radio power which 
can be interpreted as evidence for a low electron acceleration 
efficiency. In particular, we find that £ e < 0.001. According 
to the magnetic scaling proposed bv lBonafede et all l|201ol i 
in the case of the Coma cluster (model 'b') the acceleration 



efficiency could easily reach values of £ e 



10" 



However, 



there are many uncertainties which may affect the Po value, 
e.g. the actual discovery probability. 

• After normalizing the radio relic number counts, 
N(> S u ), we split the obtained number counts into redshift 
bins. As a result we expect more relics for z > 0.3 than 
being observed. This might indicate that the B(n B ) relations 
assumed in the simulation show an additional dependence 
on redshift. However, magnetic model 'b' seems to agree 
better with observational results which would point toward 
> fiG magnetic fields in the location of radio relics. We 
consider this approach as a very promising diagnostics of 
the evolution of magnetic fields in galaxy clusters but larger 
relic samples are needed to draw any robust conclusion. 

• The observed relic number counts are reasonably 
reproduced assuming an effective sensitivity of 100 mjy. 
Candidates for many relics may have been first identified in 
the NVSS survey, which has a r.m.s. noise level of 0.45 mjy. 
We adopt therefore that the effective sensitivity for finding 
relics is generally about 200 times the r.m.s. noise of a 
survey. We apply this to the specifications of the proposed 
LOFAR and APERTIF surveys. Under these assumptions 
we find that the LOFAR- Tier 1-120 MHz has the potential 
to find more than a thousand radio relics. 

• More than 50% of the relics expected to be found with 
LOFAR- Tier 1-120 MHz survey should reside in clusters 
with z > 0.3 and there should be even more than 100 relics 
in clusters with z > 1. Hence, in principle this survey will 
allow to discover sufficient relics to analyze the evolution of 



magnetic fields in clusters in a statistical way. 

• To confidently discover a relic, the clusters which host 
the diffuse radio emission need to be identified. Many of the 
relics which are detectable by the LOFAR-Tier 1-120 MHz 
survey may reside in faint clusters. More precisely, we 
predict that about 50% of the relics with 20 mjy will reside 
in clusters with an X-ray flux below 4 x 10~ 13 ergs -1 cm -2 . 



• Following van We eren et all l|2011ch we study an X-ray 
flux galaxy cluster sample based on the NORAS+REFLEX 
catalogues. About 4% of the galaxy clusters in the sample 
host radio relics. This value is significantly lower than the 
one obtained for the overall sample since we consider only 
known clusters here. We found that we can reproduce the 
relic fraction assuming a n effective sensitivi t y of 27 mjy at 
1.4 GHz. As discussed in Ivan Weeren et all feoilcft we also 
find that fraction of clusters which host a relic increases 
with cluster X-ray luminosity and redshift. 

• We expect that the LOFAR-Tier 1-120 MHz survey will 
find radio relics in around 50% of the NORAS+REFLEX 
clusters. Furthermore, for the most massive clusters this 
fraction can be as high as 90%. 
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